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The rationale for introducing non-hermitian Hamiltonians and other observables is reviewed and 
open issues identified. We present a new approach based on Moyal products to compute the metric 
for quasi-hermitian systems. This approach is not only an efficient method of computation, but also 
suggests a new perspective on quasi-hermitian quantum mechanics which invites further exploration. 
In particular, we present some first results which link the Berry connection and curvature to non- 
perturbative properties and the metric. 
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I. INTRODUCTION 

The history of quantum mechanics now spans more then a century. The mathematical underpinning, although 
more recent, is also well established, particularly due to the work of von Neumann and Weyl, amongst others (see [lj 
for a comprehensive treatment). A central result of this analysis, known as the Stone- von Neumann theorem (yj, vol 
2, pp. 34-37), states that for systems with a finite number of degrees of freedom, all unitary irreducible representation 
of the Weyl algebra (determined by the canonical commutation relations) are equivalent. In essence this implies that 
the quantization of systems with a finite number of degrees of freedom is unique, up to unitary transformations. For 
systems with an infinite number of degrees of freedom the situation is more complex. In this case there are inequivalent 
representations of the Weyl algebra and, correspondingly, such systems can be quantized in inequivalent ways. 

Central to the analysis mentioned above, which forms the basis for our current formulation and interpretation of 
quantum mechanics, is the idea that classical observables should be represented as hermitian operators on the quantum 
level. In particular position and momentum are viewed as observables that have to be quantized as hermitian operators. 
In the past three to four decades, it has, however, transpired that this may be an unnecessarily restrictive point of 
view for a variety of reasons. In particular it turns out that it is convenient to resort to non-hermitian observables, 
mostly the Hamiltonian, when considering effective interactions 0, Q or when bosonizing a fermionic system, as the 
one- plus two-body character of the Hamiltonian can only be maintained in this way 0, Q , while also maintaining a 
lowest order association between bosons and fermion pairs. 

More recently, there has been considerable interest in so-called PT-symmetric quantum mechanics in which the 
Hamiltonian is taken to be non-hermitian, but PT-invariant. Several examples of such Hamiltonians have been found 
which exhibit a PT-unbroken phase in which the eigenvalues are real and a normal quantum mechanical interpretation 
is made possible by the introduction of a new inner product, based on the so-called C-operator having properties very 
similar to the charge conjugation operator jflQ- The possibilty of quantum computing has necessitated the study of 
open quantum systems in which the Hamiltonian necessarily becomes non-hermitian in order to capture the dissipative 
nature of the system, see e.g. Therefore, the use, and even necessity, of non-hermitian Hamiltonians is nowadays 

generally accepted, be it for mathematical simplicity or the description of physical reality 

Above we alluded to two quite distinct situations in which non-hermiticity is introduced. In the case of dissipative 
systems the Hamiltonian will generally have complex eigenvalues and a normal quantum mechanical interpretation is 
simply not possible, and indeed inappropriate. In the cases of bosonization, effective interactions and PT-symmetric 
quantum mechanics, non-hermiticity is introduced with other considerations in mind, such as a simplified mathematical 
description. A normal quantum mechanical interpretation should, however, still be possible. In this article we consider 
the latter systems, and consider when and how a system of non-hermitian observables constitute a consistent quantum 
system. This question is not new and has, in fact, been considered in considerable detail some time ago by the present 
authors y~(J. The conclusion that came out of that analysis was that a system of non-hermitian observables forms a 
consistent quantum system if a new inner product exists with respect to which the observables are all hermitian. The 
existence of this new inner product was tied to the existence of a positive definite, bounded hermitian metric operator 
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O, defined on the whole Hilbert space, and which obeys the following exchange rule with the observables Ai [lOj: 

A t Q = &A\ Vi . (1) 

Furthermore it was shown that this metric operator is unique (up to a global normalization) if and only if the set of 
observables under consideration forms an irreducible set |lOj . All these results apply to the case of a finite number of 
degrees of freedom. 

The considerations above open up an alternative way of thinking about quantum mechanics. Conventionally 
one chooses a Hilbert space (and by implication inner product) and the construction of observables is dictated 
by the requirement of hcrmiticity with respect to this inner product. With the results above in place, one can, 
however, take the point of view that the observables are the primary objects which uniquely fix the Hilbert space 
and inner product. The latter, of course, only exists if the observables form a consistent set. To implement this 
approach in practice, as is demonstrated in section UTTl one starts with a Hilbert space which carries an irreducible 
unitary representation of the Weyl algebra, i.e., hermitian position and momentum operators and thus a unitary 
representation of the Heisenberg algebra. Then one proceeds to write down the observables of the theory, e.g., 
energy angular momentum etc. as functions of the Hermitian position and momentum operators, but, in contrast 
to the conventional approach, the condition that these observables are hermitian w.r.t. the inner product on the 
Hilbert space is relaxed. To verify that this set of observables constitute a consistent quantum system one verifies the 
existence of a metric operator with the properties stated above. One can also take a slightly different point of view 
to this construction. If a metric operator exists, one easily verifies from Q that under the similarity transformation 
Q^ 1 / 2 A i (x,p)Q 1 ^ 2 = A(9 _1 / 2 a;6 1 / 2 , 6~ 1 / 2 p6 1 / 2 ) the observables become Hermitian on the Hilbert space under 
consideration (see also 0,^2)- Thus, one could take the point of view that the observables are hermitian functions 
of non-hermitian position and momentum operators 8~ 1 / 2 i0 1 / 2 and 9 _1 / 2 p0 1 / 2 . From this point of view one thus 
considers a non- unitary representation of the Weyl algebra, i.e., non-hermitian position and momentum, and construct 
the observables as hermitian functions of these non-hermitian position and momentum operators. If the metric exists, 
this non-unitary representation of the Weyl algebra is equivalent to a unitary representation. This also implies that 
if the non-unitary representation under consideration is not equivalent to a unitary one, the metric cannot exist. 

Note that since all irreducible unitary representations of the Weyl algebra are equivalent, all non- unitary represen- 
tations, equivalent to a unitary representation, are also equivalent. Whether the Weyl algebra admits non-equivalent 
non-unitary irreducible representations is, to our knowledge, an open question. However, if such representations were 
to be encountered in the construction procedure above, the metric cannot exist or the construction procedure has to 
be revised. Thus, at least within the construction procedure outlined above, no incquivalcnt quantizations of a given 
system (with finite number of degrees of freedom) can occur if it is to constitute a consistent quantum system. 

Given this situation, one may wonder what is to be gained in quantum mechanics from taking the point of view 
advocated above. One advantage has already been mentioned, namely, mathematical simplification. As already 
pointed out, when bosonising fermionic systems the one-plus- two body character of the Hamiltonian can only be 
maintained in a non-hermitian description, while the hermitian counterpart contains all possible higher order inter- 
action terms 0,0. Related to this observation, it was recently also noted in the context of PT-symmetric quantum 
mechanics that the non-hermitian Hamiltonian often lends itself more readily to a perturbative treatment than its 
hermitian counterpart 13]. This reflects the fact that the implied similarity transformation relating the Hermitian 
and non-hermitain Hamiltonians already accounts for some non-perturbative effects. These examples illustrate that 
a judicious choice of observables (which fixes the inner product) may lead to considerable simplification and even the 
possibility of non-perturbative solutions. 

Central to the succesful practical implementation of the approach outlined above, is the verification of the existence 
of the metric operator and its subsequent construction. In this regard it is important to note that if one is only 
interested in solving the eigenvalue equation for the Hamiltonian, it is sufficient to verify the existence of a metric 
operator obeying the exchange rule with the Hamiltonian as one is then ensured of the reality of eigenvalues. For 
this purpose it is not necessary to construct the metric explicitly as the solution of the eigenvalue equation requires 
no explicit knowledge of the metric. However, if one intends to compute other physical quantities such as expectation 
values, transition probabilities and cross sections, detailed knowledge of the metric is required and the metric must be 
computed explicitly. It is also at this point that the uniqueness (up to a global normalization) of the metric becomes 
important as physical quantities will depend on the choice of metric if it is not unique. 

Solving operator equations of the type (JIJ is notoriously difficult. The succesful implementation of the approach 
outlined above hinges therefore very much on a reliable way of solving these types of equations, or at the very least 
verifying the existence of a solution with the properties required from a metric operator. Therefore this issue was 
partially investigated in 01 an d more recently from a variety of different viewpoints in [Til ITU [Til IT7L fTsL Il9| . Here 
we propose a new approach, based on the Moyal product, which maps the operator equation exactly onto an 
equivalent differential equation. Generically this equation may be of infinite order, but in many cases of physical 
interest it turns out to be finite. This equation contains all the information required to construct the metric operator 
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exactly. In addition, criteria can be formulated to test the hermiticity and positive definitness of the metric directly 
on the level of this equation, leading to considerable simplification. 

When one considers eq. for a reducible set of observables, e.g. the Hamiltonian alone, the metric is not uniquely 
determined. However, on the level of the Moyal product formulation the solution of the corresponding differential 
equation is uniquely determined once an appropriate set of boundary conditions has been specified. On the other 
hand, as pointed out above, the metric is uniquely determined (up to an irrelevant normalization factor) if cq. Q 
is considered for an irreducible set of observables. This suggests an interplay between boundary conditions in phase 
space and the choice of an irreducible set of physical observables. This provides us with yet another perspective on 
the construction of a quantum system, which is the one we explore somewhat more here. 

The paper is organised in the following way. In section [n] we briefly review the Moyal product with emphasis on 
the properties we use in later sections. In section UTTl we construct the metric of three non-hermitian PT symmetric 
Hamiltonians, two of which are exactly solvable. The issues of uniqueness of the metric, choice of observables and the 
relation to boundary conditions in phase space are addressed within these models. In secti on [VI we conclude with a 
brief summary of our findings. The present paper extends the results recently published in |2f| . 



II. THE MOYAL PRODUCT 



The Moyal product [2l| was introduced by Moyal to facilitate Wigner's phase space formulation of quantum me- 
chanics. More recently it was revived in the context of non-commutative systems (see e.g. H^)- In this section 
we briefly list the essential features of this construction. A more detailed exposition can be found in (2jj. In a 
finite dimensional Hilbert space with dimension N one can construct an unitary irreducible representation of the 
Hcisenberg-Weyl algebra 

gh = e i *hg; = g~ l , h) = hT 1 , (2) 

where <j> = The operators U(n, m) = g n h m , with n = 0, 1, . . . N — 1 and m = 0, 1, . . . N — 1, form a basis in the 
space of operators (matrices) on the Hilbert space. Any operator A can therefore be expanded in the form 

JV-l 

A= ]T a n , m g n h m , a n , m = (U(n,m),A)/N. (3) 

n,m— 

Making the substitutions g — > e la h — ► e l{3 , a , (3 G [0,27r) in the expansion turns A into a function A(a,f3), 
uniquely determined by the operator A 

N-l 

A= a n , m e ma e m P. (4) 

n,m— 

An isomorphism with the operator product can be established by defining the Moyal product of functions A(a, 0) 
and B(a,j3) [HH^ 

A(a,/3)*B(a,/3) = A(a, {3)e i(t > d ' 3da B(a, (3) , (5) 

where the notation d and d denotes that the derivatives act to the left and right, respectively. On this level operators 
are replaced by functions, as described by |@J, while the non-commutative nature of the operators is captured by the 
Moyal product. It is easily checked that the Moyal product is associative, as one would expect from the associativity 
of the corresponding operator product. Once the function A(a, j3) is given, the coefficients a n>m are computed from 
a simple Fourier transform. Insertion of these coefficients in © enables the reconstruction of the operator. 

Let A'(a,f3) and A(a,f3) denote the functions corresponding to the hermitian conjugate operator A^ and the 
operator A, respectively. One easily establishes the following relation between these functions [2(1 ] 

A*{a,0) = e i * d ° a ? , A*{a,/3). (6) 

From this it follows that an operator is hermitian if and only if 

A* (a, (3) = e- 1 * 9 ^ A{a, 13) . (7) 
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These results are easily generalized to the case of an infinite dimensional quantum system. In this case a well known 
irreducible unitary representation of the Heisenberg-Weyl algebra exists [ij 

e m e isi = e iUs e isx e m ; (g) 

where x and p are the hermitian position and momentum operators satisfying canonical commutation relations. The 
operators U(t : s) = e ltp e lsx constitute a complete set pj and any operator can be expanded as 

r°° . . „ ft 

A(x,p) = / dsdt ait, s)e ltp e lsx , a(t, s) = — (U(t, s), A) . (9) 
J-oc 27r 

Replacing x and p by real numbers turns A(x,p) into a function A(x,p), uniquely determined by A 

/oo 
dsdt a(t, S y tp e isx . (10) 
-oo 

An isomorphism with the operator product can be established by introducing the Moyal product of functions 

A(x,p) * B(x,p) = A{x,p)e in9 ^B{x,p) . (11) 

On this level we again work with functions, rather than operators, while the non-commutativity of the operators 
is captured by the Moyal product. As before associativity is easily verified. Once the function A(x,p) has been 
determined, the function a(t, s) is determined from a Fourier transform. Insertion into the expansion @ recovers the 
operator A{x,p). 

As before, let the functions A^(x,p) and A(x,p) denote the functions corresponding to the hermitian conjugate 
operator A j; (x,p) and the operator A(x,p), respectively. One easily establishes the relation p(il ] 

A\x,p) = e thd * a ?A*(x,p). (12) 

This implies that an operator is hermitian if and only if 

A*(x lP ) = e- m9 ^A(x lP ). (13) 

In what follows we shall often encounter situations where the operator A is a function of x, or p, only. It is therefore 
worthwhile to consider this situation briefly. Consider the case where A(p). From @ we have 

s ) = 7^(U(t, s),A(p)) = ^tr(e- is& e- itp A(p)) = S(s) J ^A(p)e~^ . (14) 

Substituting this result in llUfl we note that the function A(x,p) corresponding to the operator A(p) is just A(p), i.e., 
we just replace the momentum operator by a real number. Clearly, the same argument applies to A(x). 

An approach related to the one we discuss here was developed in [2^, although in that case the position and mo- 
mentum operators are used as a basis to expand the operators. Compared to the present approach, the unboundedness 
of the position and momentum operators complicates the proof of completeness. Secondly, the product rule of these 
operators is not as simple as that of the Weyl algebra. This complicates the implementation on the level of classical 
variables. The current approach therefore seems to be more generic. 



III. METRICS FROM MOYAL PRODUCTS 



A. Metric equation in Moyal form 



As was already pointed out in the introduction, there are a variety of reasons why the Hamiltonian, and other 
observables, of a system may be non-hermitian with respect to the inner product on the Hilbert space on which the 
system is quantized. The central question is then whether a consistent quantum mechanical interpretation remains 
possible. This was answered in [lOj where it was pointed out that a normal quantum mechanical interpretation is 
possible if a metric operator O exists which has as domain the whole Hilbert space, is hermitian, positive definite and 
bounded, and satisfies the equation 



He = e# f , 



(15) 
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where H denotes the Hamiltonian of the system. Once the existence of such an operator has been established, a new 
inner product can be defined with respect to which the Hamiltonian is hermitian and a standard quantum mechanical 
interpretation is possible. However, as was pointed out in |loj |. and also explored in 0,113], the condition H15JI is 
not sufficient to fix the metric uniquely, which implies that the quantum mechanical interpretation based on this 
metric, and the associated inner product, is ambiguous. The metric is uniquely determined (up to an irrelevant global 
normalization) if one requires hermiticity of a complete irreducible set of observables, At (of which the Hamiltonian 
may be a member), with respect to the inner product associated with O, i.e., it is required that (|15|l holds for all 
observables |Toj : 

A t Q = QA\ Vi . (16) 

From this point of view the choice of observables determines the metric and Hilbert space of the quantum system 
uniquely. Alternatively, one may argue that if a metric, satisfying 1)150. can be found, a particular choice of metric 
determines the allowed set of measurable observables. This is the spirit of PT-symmetric quantum mechanics. 

Let us consider the defining equations (|15f) and (|16|) on the level of the Moyal product formulation. On this level the 
observables A i} their hermitian conjugates A\ and the metric operator O get replaced by functions Ai(x,p), A\(x,p) 
and Q(x,p) as prescribed in (|10fl . Note that Ai(x,p) and A\{x 1 p) are related as in i|12|) . In terms of these functions 
the defining relation l|15fl then reads 

H(x,p) * Q(x,p) = Q(x,p) * rf(x,p) , (17) 

while lfR)|l reads 

Ai(x,p) * Q(x,p) = e{x,p) * A\{x,p) . (18) 

If the set of observables Ai is irreducible, equation (|18f) determines the metric uniquely up to a global normalization 
factor, without the necessity of additional boundary conditions. On the other hand, equation i|17|) does not determine 
the metric uniquely and it is necessary to impose further boundary conditions to fix the metric up to a possible 
normalization factor. This demonstrates the already mentioned interplay between the choice of observables and 
boundary conditions on the metric function 0(x,p). To demonstrate this point more clearly we consider the case 
where we specify x and p as observables. In this case equation (|18f) simply becomes O^ 1,0 ^(x,p) = O^ 0,1 ^(x,p) = 0, 
i.e., the metric is just a constant. This simply reflects the following facts: (1) that x and p form an irreducible set, 
which implies that the metric is uniquely determined up to a global normalization factor and (2) that we have chosen 
x and p hermitian from the outset (insisting on unitary representations of the Heisenberg-Weyl algebra), so that the 
metric must be proportional to the identity (the original inner product corresponds to O = 1). 

B. A solvable model 

To demonstrate the technique described above, and to illustrate the features discussed thus far, we study in this 
section a model for which the metric can be computed exactly and analytically. We proceed by first choosing the 
Hamiltonian as the only observable, leaving all other possible observables unspecified. It was already pointed out that 
this does not fix the metric uniquely and this is demonstrated explicitly below. Once this has been established, we 
proceed to augment the Hamiltonian with other choices of observables which, together with the Hamiltonian, form 
an irreducible set and show that this indeed eliminates the non-uniqueness of the metric. 

The model we study is defined by the following Hamiltonian, written in second quantized form: 

H = fiuja) a + Tiaaa + fi(3a) a) , (19) 

where a and a? are bosonic annihilation and creation operators, respectively. A finite dimensional version of this 
model was studied in [10J and more recently this model was studied in |18l Il9l |23| . This can be rewritten in terms of 
x and p in the usual way by setting cl> = (x — ip) / and a = (x + ip) /Vzh. Suppressing an irrelevant real constant 
term, which plays no role in the metric equation, this yields 

H = ap 2 + bx 2 + icpx , a = {uj - a - (3)/2, b = (tu + a + /3)/2, c = (a - 0) . (20) 

On the level of the Moyal bracket formulation this becomes 



H(x,p) — ap 2 + bx 2 + icpx; H\x,p) = ap 2 + bx 2 — icpx + ch . (21) 
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Substituting this in (|17fl and noting that since H(x,p) is polynomial the Moyal product truncates, one obtains the 
following partial differential equation for <d(x,p): 

c (K — 2 ipx) Q{x,p) 

+h ((cp-2ibx) Q^^ix^p) + {cx + 2iap)e {1 ' 0) (x,p) + bhQ^°' 2 \x,p) - ahe^'°\x,p)^ = 0, (22) 

where the notation ©( n > m ) = d n+m Q/d n xd m p has been introduced. 

Note that since no boundary conditions have been specified, the solution is not unique. On the level of the Moyal 
product the non- uniqueness of the metric therefore resides in the freedom to specify the boundary conditions in l|22|l . 
It should, however, be noted that the boundary condition that is to be imposed on (|22|l can not be arbitrary as the 
solution has to conform with the conditions of hermiticity and positive definiteness of the metric. Keeping in mind 
that the metric is uniquely specified once a complete set of observables, hermitian with respect to the inner product 
associated with 0, is identified, an interplay between the boundary conditions imposed on (1221) and the choice of 
physical observables on the operator level is suggested. In this regard, note that a choice of boundary conditions 
that does not admit a solution conforming to hermiticity and positive definitess, constitutes an inconsistent choice 
of observables, and subsequently an inconsitent quantum system, as was pointed out in |1C| . On the other hand, 
if an appropriate choice of boundary conditions, which fixes the metric uniquely, is made, both the Hilbert space 
and allowed observables in the quantum theory are uniquely determined. The allowed observables can indeed be 
computed by solving (|18|l for the observables once the metric has been determined. In this case each choice of 
boundary condition for solving (|18|l corresponds to an admissable observable. Of course, it remains to be determined, 
typically a posteriori, when a set of such observables, now by construction quasi-hermitian with respect to the same 
metric which has been uniquely fixed (or chosen) with respect to the Hamiltonian, constitutes an irreducible set 
together with the Hamiltonian. 

Before proceeding to the detailed solutions of (1221) , it is useful to comment on some general features of the equation 
and its solutions. The first issue that we address is the existence of an hermitian metric operator, O, as required 
by the definition of the metric operator. Equation l|22[) is clearly linear and of the form LQ(x,p) — 0, with L 
a differential operator. Using e- lhd *"p xe ^d a 9 p _ x _ faQ^ anc [ e -ihd x d v ^ e ihd x d p _ ^ _ if l g x ^ one easily verifies 
e -ihd x d p L e ihd x d p _ -X*^ so that L* e^ lhdxdp( d(x,p) = 0. On the other hand the complex conjugate of (|2*2)l reads 
L*Q*(x,p) = 0. Provided that the boundary conditions imposed on 1221) also satisfy (|13|) . uniqueness of the solution 
f (|22|) is linear) implies O* (x,p) = e~ thd:cdp ®(x,p). Thus, provided that the boundary conditions imposed are consistent 
with l|13|l . the solution of 122|) . when employed to construct the metric operator as described in the previous section, 
will yield an hermitian metric operator. 

A further property to note is that if Q(x,p) is a solution (not necessarily corresponding to an hermitian and positive 
definite operator) of (|22(l . or equivalently (|17|l . then for arbitrary functions f(H(x,p)) and g{H^ (x,p)) the following is 
also a solution: f(H(x,p)) * Q(x,p) * g{H^ (x,p)), where the functions / and g are defined through a Taylor expansion 
involving the Moyal product. This is most easily checked directly on the level of equation (|17|l using the associativity 
of the Moyal product and the fact that f(H) * H = H * f(H) and g{H^) * = * g{H^). This is, once again, 
merely a reflection of the non-uniqueness of the solution of l|22[) . which has to be eliminated through some appropriate 
choice of boundary conditions. As was pointed out earlier the boundary conditions cannot be completely arbitrary, 
but must conform to hermiticity and positive definiteness. This does not, however, eliminate the freedom in 1221) 
completely and more input is required in the form of boundary conditions. Indeed, one can easily verify from l|17|) . 
associativity, and the relation l|12|) . that if 0(x,p) is a solution corresponding to an hermitian and positive definite 
operator, so will be g(H(x,p)) * <d(x,p) * e thd:cdp g(H(x,p))* . 

Let us now proceed to the detailed solution of (1221) . It is not difficult to find an exact solution in the form of the 
following one parameter family of metrics: 

0(x,p) = e rp2+spx+tx \ (23) 

where 

-c ± Jc 2 - Aabhs (2i - hs) c ± Jc 2 - Aabhs (2i - hs) , . 

r - UH "' t= ^ 4^ "» (24) 

s being a free parameter. Note that the solution has an essential singularity at % = 0, so that the metric does not exist 
as a classical object. Once this solution has been found, a large class of solutions can be constructed as was pointed 
out earlier. The non-uniqueness of the solution is, however, already explicit in l|23|) and we concentrate on that for 
the moment. To see how this freedom can be eliminated through the specification of other observables, consider 
the situation where we specify the momentum as a further observable. Then, from (|18f) . Q(x,p) also has to satisfy 
the equation Q^'°\x,p) — 0, i.e., it can not depend on x. To satisfy this condition we must have s = 0, t = and 
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r = — 251 = (u+a+0)h ' whicli removes the freedom in l|23|) . Furthermore it is clear that the arbitrariness in the solution 
g(H{x,p))*Q(p)*e l " nda:dp g(H(x,p))* is also removed in that the introduction of the function g will lead to an unwanted 
position dependency of the metric through the position dependency of the Hamiltonian. Similarly if one specifies 
the position as another observable, the metric must be independent of momentum, which is the case when s = 0, 
r = and t = tSf = -, — ° l ~? mf . As before the introduction of the function q introduces an unwanted momemtum 
dependency in the metric, so that all freedom, apart from a global normalization factor, has been eliminated. Note 
that both these solutions also have the correct asymptotic behaviour in the limit c — > 0. As a final remark we point 
out that the first of these solutions would have resulted by imposing the boundary conditions that on a line p = pa 
the metric is a constant, while the second results from imposing the same condition on a line x = xq. 

Let us now consider the hcrmiticity and positive definitcness of these solutions. Since a, b and c are real, the 
condition (| 1 31) is trivially satisfied in both these cases and the metric is hermitian. To show positive definitcness one 
has to verify that the logarithm of the metric is hermitian. In the Moyal product formulation this implies that one 
has to find the function corresponding to the logarithm of the metric operator and verify that it satisfies (|13|l . i.e., 
one has to find the function r}(x, p) such that 1 + rj + -^ij * r/ + -^r] * rj * rj * + . . . = Q . In this case it is, however, obvious 
that the Moyal product reduces to an ordinary product so that the function corresponding to the logarithm of the 
metric operator is simply the logarithm of J22J), which is — cp 2 /2bh or cx 2 /2ah, depending on the choice of additional 
observables p or x. This trivially satisfies (|13|l so that the metric is positive definite, although not necessarily bounded. 

As a final example of how the non-uniqueness is removed through the specification of other observables, consider the 
case where the number operator N = a^a = ^-(p 2 + x 2 ) is specified as the other observable. On the level of the Moyal 

product N gets replaced by N = ^-(p 2 +x 2 ). From (|18fl and (|21[1 we then note that <d(x,p) has to satisfy an additional 
equation, which is identical to 11221) with the specific choice of parameters a = b = and c = 0. The solution of 
this equation is clearly also of the form i|23[l where r = t = ±i^jTis(2i — hs)/2fr. Since the metric has to satisfy 

both 1221) and this equation, the value of the free parameter s is fixed to be s± — ± ^ c \^-b\ ) = Ya^p\ ) 

and r = t = 2 h(a-b) = ~ 2h(a+p) ■ Only the s_ solution yields the correct asymptotic behaviour for Q(x,p) in the 
hermitian limit c — ► or a — ► /3. One can also easily verify that no further freedom exists by introducing the function 
g. The reason is simply that these functions can only depend on either H or N if one wants a solution of either of 
these two equations. Since both equations have to be satisfied this freedom is thus eliminated. 

Finally we check the hermiticity and positive definiteness of this metric, i.e., we check whether the metric and its 
logarithm satisfy (|13fl . We first check the metric. Since it is an exponential it is very difficult to verify the hermiticity 
directly. Instead, the way we do it is to check the hermiticity order by order in a series expansion in c. As an example 
we show the calculation to third order in c. To this order we have for Q(x,p) 

c(p 2 +x 2 ) c 2 {p 4 + Aihpx + 2p 2 x 2 + x 4 ) c 3 (p 2 + x 2 ) (p 4 + 12ihpx + 2p 2 x 2 + x 4 ) 
e{x > P > = 1+ 2h(a-b) + 8(a~b) 2 h 2 + 48 (a - bj h 3 

= 1 + cai + c 2 a 2 + c 3 a 3 . (25) 

It is now simple to verify that Ijl3(l is indeed satisfied to this order in c, verifying the hermiticity of the metric. Next 
we turn to the logarithm of the metric. To compute the logarithm to third order we simply expand the logaritm of 
(I25|l . keeping in mind that the expansion must be done through the Moyal product in order to reflect the operator 
nature correctly. One easily finds 

(logO)(x,p) = log(l + C ai+c 2 a 2 +c 3 a3 ) =a lC + c 2 (^a 2 - ^y^J +c 3 f a 3 + 1 - -1 M . (26) 

Substituting from l|25|l the values for ax, a 2 and 03, the logarithm of the metric can easily be evaluated to this order 
to be 

t^x*" - + (achy + W^w)) <f + x ' ] ■ (27) 

It is now simple to verify that H13(l is satisfied by 1)27(1 . verifying the positive definiteness of the metric. We note that 
11271 has a simple linear form in N = p 2 + x 2 . It is not too difficult to see that this is indeed a feature that persists 
to all orders so that in this model the logaritm of the metric has the simple form (logO)(x,p) = a(c) + b(c)N where 
a(c) and b(c) are real and a(0) = b(0) = 0. The condition l(13|l can then in fact be verified to all orders in c. This 
confirms the positive definitcness and hermiticity to all orders. 
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C. A shifted oscillator: the ix potential 

Not unexpectedly an exact solution can be obtained for the metric (and equivalently for the C-operator) when 
the Hamiltonian is the PT-symmetric complex shifted harmonic oscillator with H = ^p 2 + ^x 2 + ix. In this case it 
has been shown that the general form of the C-operator, C = e^V, simplifies to C = e 2p V . It is well known that 
the metric can be related to a similarity transformation that transforms a given non-hermitian Hamiltonian into and 
equivalent hermitian form (see also section^), and that the similarity transformation is also related to the C-operator) 
jlOt fl4t fl5j . which allows the identification = e _2p in the present example. 

This result also immediateley follows from the general analysis in section |H] In particular, eq. l(T7)l reduces to the 
partial differential equation 

2ixO{x,p) + {ix - 1)QM - i6 (0 ' 2) - iQ&Vp + ±© (2 ' 0) = 0, (28) 

where we have set K = 1. Assuming 0(x,p) = 0(p) only, eq. (|28|) becomes the ordinary differential equation (primes 
indicating differentiation with respect to p) 

2ixQ + {ix - 1)6' - i6" = 0, (29) 

with the solution 0(p) = e~ 2p , as before. 

As already pointed out, the non-uniqueness of the metric determined from the star product l)17|l is associated with 
different boundary conditions. For the PDE (|28|l above, another class of solutions may accordingly be investigated by 
assuming <d(x,p) — <d(x). The differential equation reduces to a Schrodingcr equation with a linear potential, with 
solutions the standard Airy functions, here with complex argument. These solutions, to what extent they conform to 
the requirements of Hermiticity and positive definiteness, and their physical implications, will be discussed elsewhere 

m. 

D. The ix 3 potential 

Next we consider the PT-symmetric model with Hamiltonian Q 

H{x,p) =f + igx 3 . (30) 

Here we consider the Hamiltonian as the only observable and leave the other possible observables unspecified. We then 
proceed to investigate the existence of a (non-unique) hermitian and positive definite metric. Since the PT-symmetry 
is unbroken for this Hamiltonian Q we expect such a metric to exist. 

On the level of the Moyal product this Hamiltonian and its hermitian conjugate get replaced by the functions (it 
is a sum of functions depending on p and x only) 

H{x, P )=p 2 + igx 3 , #t fop) = H * (x, p), (31) 

while the metric becomes a function <d(x,p) as defined in UlUfl . Note that in this case the hermitian conjugate of the 
operator gets replaced by the complex conjugate of the function H(x,p). It is simple to see that this is a generic 
feature of functions (or the sum of functions) that depend on x or p only, as there is no phase due to the exchange of 
the operators e ttp and e lsx (see {HJ)). 

As before the Moyal product in equation (|17|) truncates due to the polynomial nature of H(x,p) to yield the 
following differential equation for <d(x,p): 

2igx 3 Q{x,p)-3ghx 2 Q {0 ' 1 \x,p)-3ign 2 xQ (0 ' 2 \x,p)+gh 3 Q^' 3 \x,p)-2inpQ (1 ' \x,p) + n 2 Q {2fi) {x,p) = 0. 

(32) 

Solutions with the correct asymptotic behaviour can be constructed by resorting to a series representation in g for 
the solutions of l|32() . This also allows us to make contact with existing literature in which a series expansion was 
used |27] |. In the same spirit as the solvable model this series expansion can be tested for hermiticity and positive 
definiteness by using the criterion (|13|1 . 

We list the result to O(g) below, as the expression becomes rather elaborate at higher order, but in principle it is 



9 



quite simple to compute the result to any order desired: 

— 21 i e^T^ g ft 4 cl(p) ie^£^ftcl(p) 21 e ±± w jL g h 3 x cl(p) ie ±J T^ L gft 2 x 2 cl(p) 
{X,P) = 16? 2~p 8~p~ 5 + 8p~ 4 + 

e *r p g?ix 3 cl(p) . . 3igft 2 xc2(p) 3gftx 2 c2(p) igx 3 c2(p) gx 4 c2(p) ie ^ gftc3(p) 

1 c2 (w + its ^2 — + '—nrz 7T- + 



icA(p) 



Ap 3 Ap 4 Ap 3 2p 2 Aftp 2p 

21ie^ gh 4 cl'(p) 21 g h 3 x cl'(p) 9 i g h 2 x 2 cl'(p) g ft, x 3 cV \p) 



16p 5 8p 4 8p 3 Ap 2 

3igft 2 xc2'{p) 3gftx 2 c2'{p) igx 3 c2'{p) 9 i e 2 ^ g h 4 cl"(p) 9 g ft 3 x cl"(p) 

V + + 2p 16P 3 8p 

gft 2 x 2 cl"{p) , 3i<??i 2 xc2'» 3gftx 2 c2"{p) | ie^g^ 4 cl^(p) | 



8p 2 4p 2 4p 

gft 2 x c2^ 3 \p) 



Ap 2 2p 



0{g 2 ). (33) 



This is the most general form of the solution where cl(p), c2(p), c3(p) and cA(p) are completely arbitrary functions of 
p. These functions are, however, restricted if one requires that the expansion (1-i-SI) satisfies the hermiticity condition 
(|13f) . Indeed, one can quite easily see that c2(p) and cA(p) must at least be real. Here we do not pursue the most 
general solution, but rather focus on a particular choice of integration constants for which this expression simplifies 
considerably. Setting cl(p) — 0, c2(p) = 1 (this brings about only a global normalization), c3(p) — and cA(p) = 
one finds to O(g) 

. (3ift 2 x 3ftx 2 ix 3 x 4 \ 

e(x,p) = l + g — — ____ + __ . 34 

v ' \ Ap 4 Ap 3 2p 2 AftpJ y ' 

This result agrees with the result in |27j ] when the different normalization of the kinetic energy term, accounting for 
the factor of i, and the different convention for the definitions of the metric (see 11151 ). which brings about a complex 
conjugation, are taken into account. With this choice of boundary conditions, even the higher order term simplifies 
considerably and one easily finds to 0(g 3 ) (setting any further integration constants zero) 



,'3ift 2 

Q(x,p) =l+g 



„2 



Ap 4 Ap 3 2p 2 A ftp 
fl08ift 5 x 108^ 4 a; 2 57ift 3 x 3 21ft 2 x 4 6iftx 5 11 x 6 



9 \ p 9 p s p 7 + p 6 + p 5 8p 4 Aftp 3 ' 32ft 2 p 2 

( 29872557 i ft s x 29872557 ft 7 x 2 7676559 i ft 6 x 3 5395599 ft 5 x 4 727299 i ft A x 5 



9 [ 256p 14 256 p 13 128p 12 + 256 p 11 + 128p 
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159489 ft 3 x 6 lA679ift 2 x 7 9207 ftx s 615 ix 9 343 x 10 Six 11 x 12 \ 4 

128 p 9 64? + 256 p 7 + 128p 6 ~ 640 ftp 5 ~ QAft 2 p 4 + 38 4 ft 3 p 3 ) + 

= 1 + ga + g 2 b + g 3 c + 0(g 4 ) . (35) 

Note, as in the solvable model, that Q(x,p) is singular in the limit ft — ► 0. 

We can now proceed to check the hermiticity of by verifying that (|13fl holds for l|35|) . It is easily verified that 
this is indeed the case up to 0(g 3 ). 

Finally, we consider the positive definiteness of 0. For this we have to show that the logarithm of is also hermitian. 
As in the solvable model the logarithm of can be computed from (I35|l by using l|26|l . This yields 

,'3ift 2 x 3 ftx 2 ix 3 

(log 8) fop) 



9* 



Ap 4 Ap 3 2p 2 Aftp , 
-2745171 if? x 2745171 ?i 7 a; 2 677457 ift 6 x 3 439857 ft 5 x 4 52029 ift 4 x 5 



256p 14 256 p 13 128p 12 256 p 11 128? 1 

9375 ft 3 x 6 651ift 2 x 7 273 ftx 8 5ix 9 



128 p 9 6Ap 8 256 p 7 64p 6 320 ftp 5 



+ onrr* r , (36) 
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which can easily be checked to satisfy ff 1 31) . Note that the second order term vanishes with this choice of boundary 
conditions 27] . We have now verified the hermiticity and positive dcfinitcness of the solution, at least to third order 
in the coupling g. 



IV. THE BERRY CONNECTION AND CURVATURE 



In this section we show how the formalism developed in section [H] can be used to construct the Berry connection 
and curvature [2^, [2^, [3(| • Although we cannot verify directly from this the existence of a metric operator for a 
given Hamiltonian, we can identify points or lines in parameter space where singularities will occur in the metric, 
eliminating the possible existence of a metric in these cases. This enables us to make statements about the analytic 
properties of the eigenvalues and eigenstates, which is a non-perturbative piece of information providing information 
on the radius of convergence of a perturbative treatment. 

Let us start by setting up the equation for the Berry connection. Consider a Hamiltonian depending on a set of 
real parameters q±, qi • • • <Znj denoted H(q). Here we are interested in non-hermitian Hamiltonians and therefore the 
Hamiltonian does not have to be hermitian as we range over the parameter space. Furthermore we do not even insist 
on real eigenvalues of the Hamiltonian as we range over parameter space, in which case a metric can of course not 
exist. We now write the Hamiltonian as 

H(q)=S(q)D(q)S- 1 (q), (37) 

where D(q) is a diagonal operator. It is important to note here that the operator S(q) may be singular at certain 
points in parameter space. These are the points where the Hamiltonian is not diagonalizable (it can at best be brought 
in the Jordan form) and the eigenstates do not span the whole space, i.e., they are linearly dependent. At these points 
the metric can also not exist (if it does, one can construct an operator that diagonalizes the Hamiltonian). To proceed 
we differentiate (|37[1 with respect to qi 

d -^-S{ q) d -^S-\q) = [MqlH{ q )]. (38) 

Here we have introduced the Berry connection defined by 

MQ) = ^S-Hq) ■ (39) 
oqi 

The Berry connection simply generates the change in the eigenstates as the parameters on which the Hamiltonian 
depend are changed. Furthermore one expects singularities in S(q) to show up in the Berry connection. Thus solving 
for the Berry connection from (|39|1 one may be able to identify these singular points. Note that as the components of 
the Berry connection are operators they do not necessarily commute. 

Equation (|38l) is not very useful in its given form, as we need to know the eigenvalues in order to compute the 
Berry connection from it. To avoid this, we take the commutator of (|38|) with the Hamiltonian. It is simple to see 

that [ S (q)°§&S-Hq),H(q)} = S(q)(D(q)°§& - d -^D{q))S-\q) = and we have 

[^.^)] = [[Mq),H(q)},H(q)}. (40) 
oqt 

We can now use this equation to solve for the Berry connection. Once again this is an operator equation which is 
generally difficult to solve. On the level of the Moyal product this becomes, however, again a differential equation of 
finite order if the Hamiltonian is polynomial. 

Before attempting to solve (|4()|> in a particular model, we first have to consider some of its general features. It 
is clear that (|40|l does not have a unique solution as we can always add a solution of the homogenous equation to 
find another solution. This freedom is indeed already present in equation (|37|l as a 'gauge freedom'. To see this we 
note that if S(q) diagonalizes H(q), so will T(q) = S(q)A(q) where A(q) is an arbitrary diagonal operator. Under the 
transformation S(q) — » T(q), the Berry connection transforms to 

Uq) -> A'M = AM) + S(q)^-A- l (q)S- 1 (q) , (41) 

oqi 

illustrating the non-uniqueness of the Berry connection. Quantities like the Berry phase should, however, be unique 
and thus invariant under the transformation Q41|). Let us therefore consider this issue in some more detail. To do this 
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FIG. 1: Change in S(q) with a translation around a plaquette in the qi and qj directions. The initial value of S(q) is denoted 
by So and the final value by 5*4- Intermediate values are denoted Si, S2, and S3, repectively. 



we consider the change in S(q) as we translate around a small plaquette in the qi and qj directions as shown in figure 
^ Expanding to second order in dqi and dqj , one easily finds the following relations between the intermediate values 
(see figure of S(q) 

Si = (l + A^ + l^+A^dq^jSo, 

dA t , , 1 fdAi 
dqj 



S 2 = ( 1 + Aid qi + -^dq t dqj + — ( -77— + Aj ] dqj ) S 1 



S,= [l- A j dqj - dqj + \ [A] - ^ ) ,Ur j S 2 

S, = li-A^ + lU-^dqAss. 



(42) 



Successive application now yields the relation between S4 and So 

Si = (1 + Fijdqidqj) S , (43) 
where we have introduced the Berry curvature 

*>-%-?£ + »*W- < 44 > 

This is the infinitesimal change in the eigenstates. For finite closed loops the change can be obtained from a path 
ordered exponential of the Berry connection or a careful exponentiation of i|43[) . keeping in mind that the at 
different points need not commute. It is now a simple matter to verify that under the transformation l|41J) the Berry 
curvature is invariant, i.e., F!j — Fij. 

Before applying l|40|) to infinite dimensional quantum systems, we consider a simple 2-dimensional matrix problem. 
Consider the following Hamiltonian 

"<^>=(„+ I51 "- + r) <«> 

One can easily solve (|4(J|I to find the most general solution 

Ai(qi,q 2 ) = ( qi+iq2 ~ (n+iQ2) (i+(?i+i« 2 ) 2 ) + Vl Wl ~ i+(9i+«?2) 2 j 

\ wi yi J 

/ 2iwi _ i 1 : _ i \ 

A 2 (qi,q 2 ) = 91+492 (9i+'?2)(i+(9i+»g2) 2 ) yi i+(<ji+i<?2) 2 ( 46 ) 

\ IWl 2/2 / 
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The homogenous part has been constructed to yield a zero curvature. We note that this expression has singularities 
at the points {ql = 0, q2 = 0} and {ql = 0, q2 = ±1}. The singularity at the origin is spurious (the transformation 
S(q) = 1 at this point) and can be removed by an appropriate choice of the homogenous part. Indeed choosing 
wi = 1/2 and y% = y 2 = we have 



<?l+'g2 1 



A 1 {q 1 ,q 2 ) = -iA 2 (q 1 ,q 2 ) = I i+(<^) 2 2 i+^+mfj (47) 

The singularities at the points {qi = 0, q 2 = ±1} are, however, not removable and, not surprisingly, are the two 
exceptional points where the matrix is not diagonalizable |3l) and the transformation 5(g) becomes singular. When 
one computes the curvature it vanishes everywhere, except at these points where the curvature has to be computed 
with great care due to the singularity. This implies that moving in an infinitesimal loop around any point brings the 
eigenfunctions back to themselves; the only exception being the exceptional points where the eigenfunctions transform 
into each other [3l]]. Note that the Berry connection at different points does not commute. Due to this fact one has 
to compute a path ordered exponential in order to obtain the change of the eigenvectors as one moves around in a 
closed loop. The result (|47|l can also be obtained from (|39|) by direct diagonalization and an appropriate choice of 
'gauge'. 

To demonstrate the aforementioned points, we compute the curvature at the exceptional point {qi =0,(72 = 1}- 
We do this by computing the change of the eigenstate when we move in an infinitesimal circle around the exceptional 
point. Setting q\ — rcoscj) and qi = 1 + rsin</>, we easily compute the Berry connection in the azimuthal direction 
from (|39(l and l|47(l to be = ire 1 ^ A\. This result is now expanded to lowest order in r to obtain 

^5-i(0)l(o.i) =^ = ( J 7) ■ ^ 

The change in the eigenstate can now easily be computed from (|48|l by noting that S(2n) = (1 + ^-A,p) N S(0) in the 
limit N -> 00. The result is 5(2tt) = FS{0) with 

. = (-?). m 

We now proceed to compute the action of F on the eigenstates infinitesimally close to the exceptional point. These 
eigenstates, expanded to leading order in r, read 

( — i ± iy/2w \ 

u ± = ( o ) > (5°) 

where we have set w = re 1 ^ . It is now a simple matter to see that Fu± — u^. 

Let us now consider the quantum system discussed in section IlII Bl All the considerations of section ITT1 apply and 
on the level of the Moyal formulation we can replace equation (|40l) by 



dH{x,p, q) 

^ ; H (x,p, q) 

oqi 



[[Ai(x,p, q), H{x,p 7 q)},H(x,p, q)}* . (51) 



Here x and p are the phase space variables, q the parameters in the Hamiltonian and we have introduced the Moyal 
bracket 

[A(x,p), B(x,p)]. = A(x,p) * B{x, P ) - B(x,p) * A(x,p) . (52) 

For the system of section UlI Bl the Hamiltonian was given in (|2H and the two parameters entering the Hamiltonian 
are the parameters a and (3 (see (|19|) and (J20I))- As we are not interested in a global rescaling of the Hamiltonian, 
which only effects the homogenous part of the Berry connection, it is more convenient to change to two new variables, 
<7i and q 2 , by dividing the Hamiltonian l|21l) by a: 

tt ( \ 2 , 2 , ■ b uj + a + P c 2(a-/3) 

H{x,p)=p +qix +iq 2 px; q x = - = g 2 = - = -. (53) 

a uj — a — p a uj — a — p 

It is now easy to set up the explicit equations for the connection from l|51|l . but we do not list them here due to 
their length. We rather make an assumption for the connection of the following form 

Ai(x,p) = np 2 + Sixp + t t x 2 ;i = 1,2 , (54) 
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FIG. 2: Singular curves of the Berry connection. These curves divide the parameter space in disjoint regions, shown as shaded 
and unshaded here, and it is not possible to pass from one region to another without crossing a singularity in S(q). The circle 
denotes the expected radius of convergence of a perturbative expansion around the origin. 



and solve for the parameters r.;, Sj and U as functions of qi and q 2 - The result is 

Ai(x,p) 



ipx x 2 q 2 



h{A qi +q 2 2 ) 2h (4 9l + g 2 2 ) ' 

M '' P) = M4 ? i+^) + 2k { 4q 1 + q 2 2 ) • (55) 

The curvature can also be easily computed from Px,2 — 7^ — + [^1,^2]*- It vanishes everywhere, except at 
possible singularities. Let us therefore investigate the singularities and interpret them. There are singularities on 
the curves 4qi + q\ = uj 2 — Aafi = 0. It is well known |2^| that on these curves the Bogoliubov transformation that 
diagonalizes the Hamiltonian becomes singular, which is now simply reflected on the level of the Berry connection. 
These continuous curves divide the a-/3 parameter space into disjoint regions and it is not possible to move from the 
one to the other without crossing a singularity of S(q). This signals that perturbation theory is limited to within 
these regions and that the radius of convergence will be set by the distance to the closest point on these curves. This 
is illustrated in figure [2J 

One may also evaluate these results from the perspective of the existence of the metric 0. While S(q) is singularity 
free, except on the border between shaded and unshaded areas, O does not necessarily exist elsewhere, since positive 
definiteness and hermiticity of 9 cannot be immediately inferred. At the same time the singularity of S(q) on the 
borders indicate the existence of zero norm states, implying that on these borders the metric definitely does not exist, 
confirming in this particular case what was anticipated in general at the beginning of this section. 



V. CONCLUSIONS 



We have presented a brief overview of considerations pertinent tot non- and quasi-hermitian quantum mechanics in 
which the role of the metric is stressed. Subsequently we have shown how the Moyal product can be used to compute 
the metric for a given non-hermitian Hamiltonian from a general partial differential equation. The verification that 
the metric is hermitian and positive definite can be carried out directly on the level of the Moyal product formulation, 
without refering to the operator level. 

We have carried through this program for three Hamiltonians, all of which possess PT-symmetry. These consid- 
erations can also be applied to finite dimensional models as those studied in [lOj, by using the finite dimensional 
formulation of the Moyal product where, essentially, h — > -k. An interesting new perspective that arises from the 
present formulation is the relation between the choice of observables and boundary conditions imposed on the metric 
equation, as formulated in terms of the Moyal product, both of which give rise to a unique metric. We have also di- 
cussed the Berry connection in the context of the Moyal product and demonstrated that non-perturbative information 
can be extracted from the resulting equation. 
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